Magnetic resonance within vortex cores in the B phase of superfluid '^He. 
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We investigate a magnetic susceptibility of vortices in the B phase of multicomponent triplet 
superfluid "^He focusing on a contribution of bound fermionic states localized within vortex cores. 
Several order parameter configurations relevant to different types of quantized vortices in "^He B are 
considered. It is shown quite generally that an ac magnetic susceptibility has a sharp peak at the 
frequency corresponding to the energy of interlevel spacing in the spectrum of bound fermions. We 
suggest that measuring of a magnetic resonance within vortex cores can provide a direct probe of a 
discrete spectrum of bound vortex core excitations. 

PACS numbers: 67.30. he, 67.30. hj, 74.25.nj 
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I. INTRODUCTION. 

Probing of a quasiparticlc spectrum in Fermi superflu- 
ids containing quantized vortices has been a challenging 
problem since the pioneering work of Caroli, de Gennes 
and Matricorti. They predicted theoretically that the 
internal electronic structure of quantized vortices in su- 
perconductors should consist of low energy fermionic ex- 
citations localized within the vortex cores with character- 
istic interlevel spacing defined as A^/Ep <^ Ag, where 
Aq is the energy gap far from the vortex line and Ep is 
the Fermi energy. For conventional s-wave superconduc- 
tors with axially symmetric vortex lines the low-energy 
del ^ Aq) part of the quasiparticlc spectrum has the 
following form 



e{fj,,kz) ~ -huvfi. 



(1) 



Here /z and k^ are the projections of quasiparticlc angular 
and linear momenta onto the vortex axis and 



huj„ 



Ao 



ykp k^E, 



(2) 



where ^ is a coherence length in superconductor or su- 
perfluid. Due to the quantization of angular motion the 
projection of angular momentum ^ takes discrete values 
which arc half integer for the s wave superconductor and 
integer for the p wave superconductors and superfluids 
like ^ He (see Refl2|). Thus the spectrum of bound vor- 
tex core fermions (tlj consists of a ladder of states which 
varies from Aq at ^ = — oo to Aq as yU = oo. An example 
of the Caroli - de Gennes - Matricon spectrum for the 
vortex in p-wavc superfluid is shown in Fig[T^. 

Bound fermions determine many of the physical prop- 
erties of vortices in superconductors and Fermi superflu- 
ids like '^He. In particular, vortex dynamics is determined 
by the kinetics of vortex core quasiparticles (see for ex- 
ample review Refy and also Refs|j-[g) and low temper- 
ature thermodynamic properties are determined by the 
peculiarities of spectrum of bound fermions. Recently 
there has been much attention focused on the vortex 



core states in chiral triplet superconductors in connec- 
tion with the topologically protected zero energy states 
being a realization of self-conjugated Majorana fermions 
(see for example RefQ). Such objects have topological 
nature^ and are the key ones for the realization of topo- 
logical quantum computation which is presently generat- 
ing much interest^. 

One of the most striking demonstrations of the exis- 
tence of such bound states have been tunneling spec- 
troscopy (STM) experiments showing a peak of the lo- 
cal density of quasiparticlc states for the scanning tip 
position at the vortex center—. However the spatial and 
energy resolution of STM techniques has been insufficient 
to investigate a discrete energy levels of vortex core quasi- 
particles which were predicted by Caroli, de Gennes and 
Matricon. 

Another experimental tool to study the discreet na- 
ture of the vortex core quasiparticlc spectrum would be 
measuring of the dynamical responses of Fermi superfiu- 
ids determined by the motion of vortices. In particu- 
lar the ac conductivity of superconductors in superclean 
regime should demonstrate a peak at the frequency de- 
termined by the interlevel transitions of quasiparticles 
localized within the vortex cores^S. The same is also 
true for the drag viscosity coefficient in the B phase of 
■^Hc (sec Refsjjlj). Of particular interest are the effects 
of fermion bound states in a vortex core on the rota- 
tional dynamics of vortices with spontaneously broken 
axisymmetry-ii. It has been predicted that a resonance 
absorption of an external rf field can occur at the exter- 
nal frequency comparable with the interlevel distance of 
localized states which is similar to the cyclotron Landau 
damping. 

However experimental realization of the mechanisms 
described above is still lacking. In superconductors the 
main obstacle is a rather restricting condition for the ma- 
terial purity, which is hardly obtainable in conventional 
superconductors. In superfluid '^Hc the reason is an ex- 
tremely high frequency of vortex motion which is required 
to excite the transitions between quantized energy levels 
inside vortex cores. Thus it seems that probing of the 
Caroli- de Gennes - Matricon theory with the help of ex- 



isting to date experimental approaches is unattainable. 

(a) s (b) 




FIG. 1: (a) Sample spectrum of Caroli-de Gennes-Matricon 
states for vortex in p-wave superfluid. (b) Sketch of tlie sys- 
tem consisting of vortex line directed along the z axis under 
the action of magnetic field. The time dependent magnetic 
field H = H(i) leads to the transitions in the ladder of Caroli- 
de Gennes-Matricon states. The interlevel transitions caused 
by the time-dependent magnetic field are shown schematically 
in the panel (a) by arrows. 

Our basic idea is to employ for this purpose the Zeeman 
interaction of quasiparticles in Fermi systems with exter- 
nal magnetic field. Quasiparticle spin degree of freedom 
has been paid little attention up to date in connection 
with the properties of the quasiparticle spectra of vortex 
cores. The reason for it is that in conventional super- 
conductors with singlet pairing the spin is a completely 
independent variable and it does not affect the orbital 
motion of quasiparticles. On the other hand the situ- 
ation is completely different in multicomponent triplet 
Fermi superfluids. Indeed, in this case the Cooper pair 
wave function is a superposition of different spin states: 



A = —i{a ■ d)(Ty, 



(3) 



where d is a vector in 3D space and (Tx.y.z are Pauli matri- 
ces in conventional spin space. In case if the direction of 
spin vector d is different at the different points of Fermi 
sphere, quasiparticles "see" the direction of Cooper pair 
spin which essentially depends on the direction of quasi- 
particle wave vector k. If there is no other spin depen- 
dent potentials then quasiparticle spin is determined only 
by the order parameter and it occurs to be dependent 
on the direction of quasiparticle propagation. Thus we 
may conclude that in multicomponent triplet superflu- 
ids quasiparticle spin and orbital degrees of freedom are 
effectively coupled. 

The schematic configuration of the system that we pro- 
pose to study is shown in FiglT}]. In general we consider 
vortices in the B phase of superfluid '^He under the action 
of external magnetic field. We assume that the magnetic 
field applied to the vortex consists of a large constant 
component along the vortex axis H^ and a small time- 
varying component directed perpendicular to the vortex 



axis H±^ = H±{t). Is we will show below due to the effec- 
tive spin - orbital coupling for quasiparticles in ■^He B the 
time dependent magnetic field can excite the transitions 
in a ladder of Caroli-de Gennes- Matrieon states through 
the Zeeman interaction with quasiparticle spins. These 
transitions are shown schematically by arrows in FigU^. 
In turn, such transitions lead to the resonant energy ab- 
sorption at the frequency determined by the interlevel 
spacing which observation would be the demonstration 
of a discrete structure of spectrum of bound vortex core 
states ^. 

It is natural to expect that the interlevel energy spac- 
ing which in '^He is of the order ujy ^ 0.1 MHz (see 
Ref. U) would determine the resonant frequency of mag- 
netic response of vortex cores. Note that such resonance 
is qualitatively different from the magnetic resonance of 
free nuclear spins (NMR) when the resonant frequency is 
determined by the constant component of applied mag- 
netic field. Within the frequency domain and range of 
magnetic fields we are interested in such type of " usual" 
NMR can be neglected. Indeed, the gyromagnetic ratio 
of "^He nucleus is of the order jg ~ 10'* Hz/G therefore 
the Larmour frequency of ul = ^gH ~' 0.1 MHz corre- 
sponds to magnetic field of about H = ^h/lg ^ 10 G 
which is much less than typical magnetic fields of the or- 
der H ^ 100 G used in NMR experiments in superfluid 
■^He. Under such conditions the frequencies that we will 
be interested in lie outside the region of usual NMR peak 
Lo^ <C ijOl which means that we can completely neglect the 
Larmour spin precession. 

This paper is organized as follows. In See.|ll]we give an 
overview of the theoretical framework which is employed 
in this work, namely the Bogoliubov-de Gennes theory 
which we use to analyze the spectra of bound fermions 
and kinetic theory to calculate a nonequilibrium mag- 
netization of vortex cores in time dependent magnetic 
field. The main results are presented in Sec. IIII[ in par- 
ticular the transformation of the quasiparticle spectra is 
discussed in Sees. IIII Al and IIIIBI The nonequlibrium 
magnetization and paramagnetic susceptibility of vortex 
cores and dissipation losses are addressed in Sec. IIII CI 
We summarize our results in Section IIVI Some of the 
details of our calculations are given in appendices. 



II. BASIC EQUATIONS. 

A. Order parameter. 

In general order parameter of triplet superfluid has the 
form given by Eq.(j31). Further assuming a, p - wave pair- 
ing which is the most relevant ease for the superfluid '^He 
we write the order parameter in the form 



Aai[l(Ja:(Jy)\ki/ k 



F ■ 



(4) 



The 3x3 matrix with complex coefficients Aai in the 
expression above can be represented as an expansion 



and finally 



A„ = A^a^^A(^A,r, 



(5) 



where A^^ = {■Xi^a±iyi,a), >H,a = ^i,a are the eigenfunc- 
tions of orbital momentum L^ and spin Sz with eigenval- 
ues /J, and V. 

Further in this paper we will deal with axially sym- 
metric quantized vortices in superfluid '^He. The axial 
symmetry of such objects is described by a generator of 
rotations around the vortex axis Q = Jz — MI, where 
Jz ~ Lz + RSz is the projection of the internal angular 
momentum of Cooper pairs onto the vortex. We intro- 
duce here the operator of 3D rotation R which trans- 
forms the coordinate axes in spin space into the ones in 
orbital space. Such rotation of coordinate axes leads to 
the transformation of the order parameter according to 
the following rule: 



^ai — RapA 



P^fli- 



(6) 



As we will see further the relative rotation of spin and 
orbital coordinate axes does not lead to qualitatively 
new results. Therefore basically we will assume that the 
spin and orbital quantization axes coincide with the vor- 
tex axis z and will just briefly discuss the appropriate 
changes in the resulting formula which take into account 
the rotation ([5]). 

The axial symmetry of the order parameter compo- 
nents satisfying the equation QAai = is described by 
the following choice of the coefficients in the expansion 

m 



C^.(r)e*(^''-^-'')'^, 



(7) 



where if, r are polar coordinates with the origin at the 
vortex center and the dependencies of the order parame- 
ter components on the radial distance r are given by the 
functions C^u{r). 

The coefficient M in Eq.(I7]) is determined by the vor- 
ticity My and the internal angular momentum of Cooper 
pairs far from the vortex axis. For example the B phase 
of '^He is characterized by the zero internal angular mo- 
mentum and therefore M — My . However in the A phase 
of ^He as well as in the chiral state of triplet superconduc- 
tor Sr2Ru04 in the homogeneous state, we have Jz = ±1 
and therefore M = M„ ± 1. 

Let us now consider the transformation of the vortex 
order parameter under the action of several discrete sym- 
metries: time inversion T, space inversion P and rotation 
by the angle tt over the x axis U2 ■ Under the time inver- 
sion we get TAai — A* ,- therefore 



Under the spatial inversion we get 



(8) 



(9) 



f/2(C^,.) = {-ir+''C-^, 



(10) 



As we will see later it is very important that we can 
construct a general expression for the pseudoscalar from 
the amplitudes of the order parameter components C^i, . 
Indeed, let us denote 



L^Ti — / ^ ^iiv 



/i+i>|— n 



Note that n = 0,1,2 since |^| = \v\ = I in ^He. Then 
a general expression of the pseudoscalar which can be 
composed from the order parameter components has the 
following gauge invariant form: 



ap = /m(Co*Ci) + /m(C2*Ci). 
From Eas. (|8l9ll0p it is obvious that 



(11) 



Pi.Cn) = (-1) 



n+M 



'-^n 



U^iCn) = i-irCn. 

Therefore assuming that the vortices are singly quantized 
M = ±1 from the above expressions we get that 

P,T,U2{ap) = -Up. 

In this paper we will consider only the vortices in the 
B phase of superfluid '^He. For singly quantized vortices 
M = ±11 there can exist five basic components of the 
order parameter. Among them are Ci__i, C_i.i and Cqo 
which correspond to the main B phase, Cq.i = Ca and 
Cifi = C/3 which correspond to the additional A and 
/3 phases localized inside vortex core. The additional A 
phase has a zero spin projection (/i = 0) and unit pro- 
jection of orbital momentum [v = 1) on the z axis while 
/3 phase has /i = 1 and v = 0^^. Far from the vortex 
core at r ^ ^„ only B superfluid phase exists so that 
(7i,_i = C*-!,! = C'oo = 1 and Ca,^ = 0. The vortex 
type is determined by the behaviour of amplitudes C^v 
at smaller distances r ^ ^y and there exist five types of 
vorticesi^. 

Vortices of o and u types are singular so that only 
the superfluid components of B phase Ci._i, C-i^ and 
Coo arc nonzero. These amplitudes are real for the most 
symmetric o vortex which is invariant under the action 
of three basic discrete symmetry transformations: Pi = 
P, P2 = PTU2 and P3 = TU2. The less symmetric u 
vortex with conserved parity Pi — P but broken P3 = 
TU2 discrete symmetry is characterized by the complex 
amplitudes of the order parameter components. 

Nonsingular v, w and uvw vortices have superfluid 
cores with the inclusion of A and /3 phases. The func- 
tions Ca,p = CA,pir) describing the spatial distributions 



of additional A and /? phases inside vortex core are finite 
at r = and vanish outside the core at r 2> ^„. The 
i; and if vortices are characterized by real B phase am- 
plitudes. If Ca,i3 are also real then we have a v vortex 
with conserved P2 = PTU2 symmetry. The case when 
Re[CA,i3) = 0, Im{CA.i3) 7^ corresponds to w vortex 
with conserved P3 = TU2 symmetry. The less symmetric 
uvw vortex with all discrete symmetries Pi, P2, P3 broken 
has complex amplitudes of B, A and /3 phases. 

Returning to the definition of the pseudoscalar pT|) 
and applying it to the particular case of vortices in '^He 
B it is easy to see that ap = for the singular o and 
u vortices as well as for the nonsingular v vortex. At 
the same time for w and uvw vortices we get a non- 
zero pseudoscalar ap ^ 0. Assuming for example the 
model situation when Ci__i = C_i^i = Cqo = CBir) 
and Ca y^ we obtain from Eq. ljlip that for w and uvw 
vortices 



ttp = Im{CgCA)- 



(12) 



As we will see below such classification of vortices in 
terms of the pseudoscalar has a direct connection with 
qualitatively different modification of quasiparticlc spec- 
trum by an external magnetic field. 



B. Quasiclassical Bogoulubov-de Gennes equations. 

Let us now turn to the spectrum of quasiparticles, 
which is described by the Bogoulubov- de Gennes (BdG) 
equations. The quasiclassical form of BdG equations 
then reads as follows: 



.hk± d 
m OS 



Ai,V= ie-P]U, 



.hk± d 
m OS 



A+C/: 



P*] V. 



(13) 



(14) 



Here k± = y^kp — fc^ is the component of Fermi momen- 
tum perpendicular to vortex axis z, P = /iB(H • a) is 
the Zeeman term and A is a gap operator ^. We as- 
sume that the magnetic field H applied to the system 
enters the equations only through the Zeeman terms and 
in general the magnetic field dependence of gap function 
Ak = Ak(H). 

Within quasiclassical formalism we should express the 
real space coordinates through the coordinate 



s = (nk ■ I") 



(15) 



along the trajectory characterized by the direction of 
quasiparticlc momentum nk = kj^/fc^ = {cos0p,sui9p) 
and the impact parameter 



(z • [nk X r]). 



(16) 



The impact parameter is related to the angular momen- 
tum projections along the vortex axis fi through the usual 
classical mechanics formula fi = —k±b. 



The energy spectrum of the quasiclassical BdG equa- 
tions £ = e(//, Op, kz) depends on trajectory angle 9p, an- 
gular momentum /i and linear momentum kz projections. 
Later we will use the following general symmetry of qua- 
siclassical spectrum (see Appendix A) 



e(/.t, Op, kz) = -e{-n, Op + tt, -kz), 



(17) 



which holds for an arbitrary triplet order parameter. 
Also we will assume for simplicity that the order pa- 
rameter is symmetric with respect to the sign change of 
magnetic field 

Afc(H) = Afc(-H) 

which is justified as long as we neglect the spontaneous 
magnetic moment of vortex core^^. Under such assump- 
tion the following symmetry of the spectrum is valid (see 
Appendix A): 



£(H) = e(-H) 



(18) 



Let us now consider the Green's functions of BdG sys- 
tem of equations (|13ll4p 






e — En as 



(19) 



where we introduce the two component eigenfunctions 
^n = {Urn Ki) of BdG equations corresponding to the 
energy level Sn- Further we will use a relation between 
the functions (|19p and the quasiclassical Green's func- 
tions g^'-^^ 

.g^(^)(s) = -z 



gR{A) 

h^kj_ 



(k, r, e) having the following form: 

G«(^)(s,s + 0) + G^(^)(s,s-0)' 



where the coordinate s along trajectory is related to the 
vectors r and k according to Ea. ([T5)) . 

Then taking the derivatives of the both sides of 
BdG system and using the normalization condition 
/ *„*+ds = 1 it is easy to obtain the expression which 
we will use later: 



/-i« 



^R ^A 



g )ds = 47r 



ft2fc 



-L V^ 



,den 



A^H^-^n)g^. (20) 



In the equation above we use the following magnetic field 
dependent part of the hamiltonian 



H 



Hb{o- ■ H) 



A^ 



Ak 

-^b(o-* • H) 



(21) 



where the diagonal terms define the interaction of nuclear 
spins with the external magnetic field and fis is a nuclear 
magneton. This hamiltonian can be rewritten with the 
help of spin operator introduced above iJ = ^b(S • H) -|- 
Ak, where we have introduced a matrix gap function 



Ak = 



Ak 

A+ 



and the operator of quasiparticlc spin 
S = {T3ax,<yy,T3az). 



(22) 



Here the Pauh matrices (^x,y,z and ti 2,3 act in spin and 
Nambu spaces correspondingly. Note that the Eq. ([20)) 
takes into account not only the energy shift due to the 
interaction of nuclear spin with magnetic field but also 
the field dependence of gap function Ak = Ak(H). In the 
next subsection we will use the relation ((20| to derive the 
expression for magnetization of vortex cores. 



C. Magnetization of vortex cores and energy 
dissipation. 

At first we need to derive an expression for the param- 
agnetic response of vortex cores in superfluid '^Hc B. Our 
essential interest is in the dynamics of vortex core magne- 
tization driven by the time-dependent external magnetic 
field and the corresponding energy dissipation due to the 
interaction of quasiparticles with the heat bath. 

We start with the exact expression of spin magnetic 
moment valid for the general non-equilibrium system: 



M 



-^J-Bi^a 



de df^k 

4 47r 



dr Tr Sg 



K 



(23) 



where we introduce the quasiclassical Kcldysh function 



g (k, r, i, e) and the operator of quasiparticlc spin 



is given by (|22[) . The integration in momentum space 
is done over the Fermi sphere so that ((ir2k)/(47r) = 
[dkzddp) / [Ai^kp) and vq is a density of states at the Fermi 
level. 

In general the average rate of dissipation losses in the 
system can be calculated as the work of the source of 
external magnetic field as follows 



ge--(M.H)t, 



(24) 



where the brackets denote the time average. However 
some parts of the total magnetization ([25)1 of the system 
correspond to a reversible exchange of the energy between 
the source and the system and for the time periodical 
processes drop out from the expression for the work of 
external source. 

To take into account such reversible energy exchange 
let us note at first that the average rate of energy loss 
does not change if we add to Eq. p4)) the time averaged 
change of the system energy, or to be more precise the 
thermodynamic potential Vl. Indeed, for the time peri- 
odical processes the average change of system thermo- 
dynamic potential is zero. Thus we can substitute the 
expression for the dissipation losses (P^ by the following 



,5Et. 



(5H 



H) 



(25) 



In the expression above we introduce the total energy 
which describes the interaction between the vortex and 



the source of external magnetic field Etot = ^ + Em: 
where D. is the thermodynamic potential of the system 
in the magnetic field and Em is the part of magnetic 
energy which describes the work done by the external 
source 5E,n = —6R = M • (5H. The variation of Etot with 
respect to the magnetic field 



5Et. 



(5H 



5VL 



M 



(26) 



gives exactly the flow of the energy to the heat bath, i.e. 
the dissipation losses. 

To evaluate the expression for the energy dissipation 
(PSj) we should calculate the variation of thermodynamic 
potential with the magnetic field 5Vt/5H. For this pur- 
pose we use the standard expression for the variation of 
the thermodynamic potential [see for example Refll4l|: 



m 



1^0 



I Tr[SHg^^^HKr)]^^dr 



^^ J(Ak(r)A*(r)) rf»k^^ 



|A| 



Att 



(27) 



where A is the weak coupling constant and the variation 
of hamiltonian ()21l) is given by 



SH 



MsS 



dA^ 



on 



sn. 



(28) 



In Eq.(|27l) we have introduced the stationary part g^'^^^ of 
the total Green's function g^ which is time independent 
and can be expressed through the stationary retarded 
and advanced Green functions g^ (k, r) and the equi- 



librium distribution function /'"'(e) 
lows: 



tanh(e/2T) as fol- 



g 



{st) 



(5, 



R 



9t)f^'\e). 



When calculating the variation of thermodynamic po- 
tential we should take into account the self-consistent 
change of the gap function determined by the following 
equation: 



Ak(r,t) = Xvq 



de driki 

4 47r 



l^(k,ki)/^(ki,r,t), (29) 



where V^(k, ki) is a pairing interaction which for the 
p-wave superfiuid has the form V^(k, ki) = (k • ki). 
Note that it is a total non-equilibrium Green's func- 
tion /^(k, r, t) which enters the self consistency equation 



Substituting the Eq. ((29| into the expression for the 
variation of thermodynamic potential (j27p we immedi- 
ately obtain that 



JO 



vq 



= f^Bl^O 



Tr [S/^*)(k,r)]-^dr 



4 47r 



Tr 



9Ak 

sn 



g("^*)(k,r) 



de df^k 

4 47r 



dr. 



(30) 



where we have introduced a non-stationary part of 
Green's function as 



9 



(nst) 



3^ 



5<"^ 



(31) 



The next step is to get use of the Eas. p3l30p substitut- 
ing them to the Ea. (|26p . Then we obtain the expression 
for magnetization ()23p to calculate the variation of in- 
teraction energy and obtain that {SEtot)/{^'H.) = Mgp 
where 



M„ 



-1^0 



4 An dU^ 



(32) 



The introduced quantity 'M.qp has the physical meaning 
of magnetization of the ensemble of quasiparticles which 
reside within vortex core. Indeed for the magnetization 
of an arbitrary system of non-interacting particles de- 
scribed by a single - particle hamiltonian H the operator 
of magnetic moment has the fornJ^ Mgp = —{dH)/{d'H.) 
and the expression p2p yields the non-stationary part of 
the thermodynamic average of 'M.qp. 

Further we will deal with the linear magnetic response 
and monochromatic processes so that it is convenient to 
introduce a paramagnetic susceptibility of quasiparticles 
in the frequency domain Xqp = Xgpi'^rf) as follows 



Mqp^UJrf) = XgpH(Wr/), 



where 



H 



Mqp{uJrf) 



/CXD 
-OO 



Mqp{t)e"^'-f*dt. 



Then for the monochromatic magnetic field the dissipa- 
tion losses are determined by the standard expression 



Qe 



Urf 



h*x;;h)„ 



(33) 



where Xqp = ImXqp- 

The expression ([5^ has an essential advantage since it 
allows to take into account the dependence of the energy 
gap on the magnetic field A^ = Ak(H), which in general 
can not be neglected. As we will see below the only thing 
we should know to calculate M^p according to Eg. ([5^ 
is a dependence of quasiparticle energy on the magnetic 
field. The quasiparticle spectrum can be either calculated 
exactly for model situations with non-self consistent gap 
Ak or it can be taken in some general form determined 
by the symmetry of the system. We will follow the latter 
way since it allows to take into account the magnetic 
field dependence of the gap function without extensive 
self-consistent calculations. 



D. Distribution function and kinetic equation. 

To proceed further we assume that the deviations from 
equilibrium are small and the rate of magnetic field vari- 
ation is much slower than the relaxation time of the gap 
function. Then we use the approximate expression for 
the non-stationary part of the Keldysh function valid for 
the slow variation of the system parametersi^: 



ff^"^*) = (5^-3^)/i 



(34) 



where the spectral Green functions g^^^^ are taken for 
the stationary system and the function /i = /i(k, r, t) de- 
termines the non-equilibrium deviation of the symmetric 
part of generalized distribution function. 

Our goal now is to rewrite the Eq. ([32|) in terms 
of the spectrum of quasiclassical BdG equation. At 
first let us use the new coordinate system defined 
by the relations (|15ll6p . Then the integration over 
drdfli^/{4:TT) in Eg. ([5^ transforms to the integration over 
dOp dfi ds dkz/{4:TrkFk±). Substituting the expression 
(pt]) into Eq.([32]) and using the relation ([20]) we finally 
obtain the following equation for the non-stationary mag- 
netization 



M 



qp 



87r2^ 



den 
d^dOpdk.—^fi 



(35) 



Finally to calculate the magnetization with the help of 
Eq. psp we only need to know the quasiparticle spectrum 
£n = ^nip-, dp) which can be found solving the BdG equa- 
tions (|13I14|) and the distribution function f{^,9p,t) = 
/o(£n(/^i ^p)) + /i: which obeys the kinetic equationii: 



5/ 
dt 



df A , df 



d9„ 



+ t/ = Stif). 



(36) 



The collision integral in the right hand size of Eg. ([55)1 
can be taken in the model relaxation time approximation: 
St{f) = if — fo)/T-i where /o is a equilibrium function 
which we take in the form: /o = tanh(e/2T). 

Canonical variables fi, Op fulfill the Hamilton equa- 
tions: 



den _ 
df-L ' 



dsn 



A* = 



(37) 



With the help of Hamilton equations the kinetic equation 
can be rewritten in the following 



(38) 



^ + {e„J}^St{f), 

where we use the Poisson bracket operator: 
df dsn df dsn 



{£«,./} = 



d9p dpi d^ dOp 



Using the Hamilton equations it is easy to see that 
{£„, /o} = as well as 

dfo \ „ 



Thus if / = /o the only term that survives in kinetic 
equation is 

dfo _ dfo dsn 
dt de dt 

Hence for the first order correction to the distribution 
function f = fo + /i we obtain the equation: 



singular vortices 



nonsingular vortices 



dt 



- + {£nJl} 



T 



dfo dSn 
de dt 



(39) 



This kinetic equation together with the expression for the 
dissipation losses (|33|) are the basic equations which we 
will use to analyze the paramagnetic response of vortex 
cores. 



III. RESULTS 



A. Quasiparticle spectrum. 



In zero magnetic field the BdG equations ()13I14|) per- 
tain the axial symmetry of wave functions which is de- 
scribed by the symmetry generator: 



Qqp = Jz~ (M/2)/T3, 



(40) 



where t^ is a Nambu spin. Then the quantum num- 
ber which characterizes quasiparticle spectrum and en- 
ters the Caroli - de Gennes Matricon expression is the 
eigenvalue of the symmetry generator Qqp"^ = /i^. 

In general due to the removed spin degeneracy the 
spectra of all singly quantized vortices in ■^He B consist 
of two different anomalous branches crossing the Fermi 
leveli^. The examples of anomalous branches for the sin- 
gular and nonsingular vortices obtained by numerical so- 
lution of BdG equations (|13ll4p are shown in Figs. ((5^) 
and ^jp) correspondingly. 

For small energies near the Fermi level an analytical 
treatment of energy spectrum of Eas. (|13ll4p is possibloi^. 
For singular o and u vortices the anomalous branches are 
similar to the standard Caroli-de Gennes -Matricon ones 
and intersect the Fermi level at zero angular momentum 
yet with different slopes corresponding to different spin 
states: 



^(m,x) = -fi^xM' 



(41) 



huj-^ ^ Ao/(fci?^) and x = ±1 corresponds to the differ- 
ent spin states. The difference in slopes of anomalous 
branches is determined by the asymmetry of amplitudes 
Ci^_i, C-i^i and Co,o inside vortex core. Further we will 
assume that this asymmetry is small to neglect the spin 
degeneracy of anomalous branches for singular vortices 
and put wi = w_i = ujy. 

On the contrary the spectral branches of nonsingular 
vortices v, w and uvw intersect the Fermi level at finite 
angular momenta [see Fig.©]: 
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H/(kp^) 





M/(kpg 



£(m,9z,x) = -fiwxAi + Ti?^ +X72, 



(42) 



FIG. 2: Anomalous energy branches e — e(/i) for (a) singular 
o and u vortices; (b) nonsingular vortices. 



where 71^2 as well as uj^ are even functions of qz = kz/kp 
and /i. Such requirements provide spectrum symmetry 
s{lJi-,qz,x) ~ ^£{^^J'^^<lz^—x) corresponding to the gen- 
eral invariance of BdG equations (jl7|) . For the nonsin- 
gular vortices it is essential that the last term in Eq. (|42)) 
corresponding to the spin splitting of energy branches can 
be rather large 72 ~ Aq. Contrary to the case of singular 
vortices when calculating the spectrum transformation 
due to the applied magnetic field we will assume that 
these energy branches are weakly interacting, i.e. the 
splitting is much larger than the Zeeman terms in Eqs. 
(IT3IT4)) . 

Now let us suppose that there is a magnetic field ap- 
plied to the system. In case when the magnetic field has a 
component H^ which is perpendicular to the vortex axis 
z the Zeeman term does not commute with the operator 
Qqp given by (j40| and its eigenvalue n is no more a good 
quantum number. Therefore the quasiclassical spectrum 
of BdG equations (|13ll4p should depend not only on fi 
but also on the conjugated angle variable which in our 
case coincides with 9p such that [Qqp, Op] = i. Note that 
the angular dependence of the spectrum is determined 
by the component of the magnetic field H^ perpendicu- 
lar to the vortex axis since the Hz component does not 
destroy the axial symmetry. 

In general the BdG system (J13I14[) can be solved only 
numerically. However the general properties of spectrum 
transformation in external magnetic field can be derived 
from the symmetry properties of the different types of 
vortices. We will consider the expansion of energy spec- 
trum by powers of magnetic field H assuming that the 
Zeeman terms are much smaller than the spacing of qua- 
siclassical levels determined by the energy scale A. 

The spectrum should obey the symmetry relations 
(|17ll8p and be invariant under the simultaneous rotation 
of coordinate axes and magnetic field around the vortex 
axis z. Then up to the first order in H± we can write the 
spectrum perturbed by the magnetic field: 

e{^i,9p,x) =eo + ai(q-H^) + a2(z- [q x H^]), (43) 

where q — 'k./kp. The first term in this expression de- 
scribes the axially symmetric part of the spectrum. Al- 
though it can depend on magnetic field for our further 
consideration this dependence will be of no importance. 



Thus wc assume that the first term in Eq.(j43]) corre- 
sponds to the spectrum without external magnetic field 
given by Eqs. ()41|42[) . The remaining terms correspond to 
the spectrum perturbation due to the field introduction 
and the coefficients ai^2 do not depend on the angle 9p. 
It is easy to check that this expression is invariant under 
the simultaneous rotation of coordinate axes and mag- 
netic field around the vortex axis z. The coefficients ai.2 
should be expressed through the characteristics of the or- 
der parameter distribution inside vortex core as well as 
the characteristics of the unperturbed quasiparticle wave 
function. 

We start with the spectrum of singular vortices which 
in zero magnetic field is given by the Eq. (|41]) with cji = 
LO-i = Uv We assume that this spectrum is degenerate 
by spin, therefore the transformation should be linear 
in Zceman shift which lifts the spin degeneracy. The 
only possibility to fulfill the requirements above is to put 
a, = xPsii where Psi are scalar coefficients, which do not 
depend on %, Op, and are even in jjl and qz- 

Now let us turn to the case of nonsingular vortices 
with lifted spin degeneracy even in non-perturbed spec- 
trum (|42p . It occurs that wc can determine the general 
properties of the spectrum (|43[) from the symmetry re- 
quirements. Indeed, the spectrum should be invariant 
under the time T and spatial P inversions and also un- 
der the rotation of the coordinate system U2- Also we 
note that the energy perturbation should be quadratic 
in magnetic field due to the spin splitting of the spec- 
tral branches (|42p and the general spectrum symmetry 
P^ with respect to the total sign of magnetic field H. 
Thus both the coefficients ai,2 should be proportional to 
the projection of the magnetic field on the vortex axis 
i^. = (H • z). 

At first let us assume that the coefficients ai_2 do not 
depend on spin. Then the only possibility to satisfy the 
symmetry requirements above is to choose 

ai,2 -- apHz (44) 

where Up is a pseudoscalar given by the Ea.(|lip: 



P,T,U2{ap 



(45) 



If we neglect the asymmetry of B phase components 
within the cores of nonsingular vortices and put C_i^i = 
Ci,_i = Co,o then from the definition of the pseudscalar 
(|12l) we immediately obtain that the spectrum of the w 
and uvw vortices then the symmetry allows that both 
the ai,2 7^ 0. 

As for the v vortex the pseudoscalar defined by Eq. 
([T^ is obviously zero in this case. However ai_2 can be 
nonzero even for the v vortex if we assume that they can 
depend of quasiparticle spin quantum number x- Indeed, 
the coefficient 72 in the spectrum (|^^ can be modified 
to include the magnetic field dependent terms as follows: 

72 = 72 (1 + Piq.HM ■ Hx) + hq.Hzi^ ■ [H^ x q])) , 

where j3i^2 are ordinary scalars. It is easy to check 
that the spectrum (|^^ with modified 72 satisfies the 



same discrete symmetries as the non-perturbed spec- 
trum. Thus the coefficients in expansion (j43p can be 
chosen as ai = xi^PiqzHz and ^2 = Xl2P2qz- 

Thus we can conclude that the quasiparticle spectrum 
of singular and nonsingular vortices has the same form 
of Eq. P5)) with a 1,2 7^ for all types of vortices. How- 
ever, the magnitude of coefficients ai^2 7^ is different 
for singular and nonsingular vortices. Indeed, in the for- 
mer case the initial spectrum is assumed spin degener- 
ate, therefore energy perturbation is of the first order in 
magnetic field, thus ai^2 ~ Ms- On the other hand the 
spin degeneracy of the spectrum of nonsingular vortices 
is lifted even at zero magnetic field. Therefore the energy 
perturbation is nonzero only in the second order of mag- 
netic field and the coefficients in (|43| are proportional 
to ai^2 ~ {lJ-%/^s)Hz, where A^ denotes the energy of 
initial split of energy branches. In general, the splitting 
can be rather large A^ ^ Ag (see Refflq) therefore we 
obtain that for nonsingular vortices the energy pertur- 
bation is smaller by the factor fisHz/^o than for the 
singular vortices. 

In general the spectrum (|43p contains the terms which 
break angular symmetry and depend on the external 
magnetic field. Hence by changing the magnetic field 
H it is possible to excite the magnetic dipole transi- 
tions of quasiparticles between the neighboring Caroli 
- de Gennes Matricon levels, which should lead to the 
resonant energy absorption for a definite frequency of 
magnetic field oscillations. 

Finally in this section we should note that the possible 
rotation of spin quantization axes with respect to the 
orbital ones given by the Eq. ^ can be easily taken 
into account in the above argument. The only thing we 
need is to transform the quasiparticle wave functions as 
follows: 



([/,F)=exp(i(a.n)^)(C/,l/), 



(46) 



where n is a rotation axis and tp is a rotation angle which 
parameterize the rotation matrix: 

{R)ik = Sik + {niUk - Sik){l- cos ip) - tikini sinip. 

It is straightforward to check that such transformation 
(|46| if applied to Eqs. (|13ll4p makes the spin axes of the 
order parameter coincide with the orbital ones. But si- 
multaneously it leads to the effective rotation of magnetic 
field 



H = i?H. 



(47) 



Note that the matrix R does not depend on the angle 
Op . Therefore the only change that should be done in the 
above consideration to take into account the rotation of 
spin axes is to replace everywhere the magnetic field by 
the rotated one (|T7)) . 



B. Numerical solution of BdG equations. 

To confirm the general argument above we solve nu- 
merically the set of quasiclassical BdG equations (|13ll4p 
to obtain the spectrum e = e{ji, Op). 

We consider the model form of the vortex core such 
that the components corresponding to the B phase are 
equal Ci._i = C_i.i = Cq.o = C's and only an additional 
A phase component is present inside vortex core. Then 
the singular part of gap operator in (s, Op) representation 
is 



Vs2 + &2 



(48) 



where 



D 



-QzCTxe 



q^a^e 



'i(l-<T-)6 



q = kp/fci? and q±^ = k^/kp^ q^ ~ kz/kp. The nonsin- 
gular part of gap function A^ is is given by 



Aa = -C^gj 



e'-^CT,. 



(49) 



In Fig. ([3]) we show the isoenergetic lines on the plane 
/i, Op corresponding to the zero energy s ~ for several 
generic cases: (i) singular o and u vortices; (ii) nonsin- 
gular vortices for H^ ~ 0; (iii) nonsingular v vortex for 
Hz ^ 0; (iv) w and uvw vortex for H^ ^ 0. 

To understand the numerical results shown in Fig. ([3]) 
let us consider the expressions for isoenergetic lines ^ = 
fi{Op, e) which can be derived from the general expression 
for quasiparticle spectrum (|43p . For singular vortices we 
obtain 

KOp) = X^i cos{Op - Oh) + xh sin(0p - 0,,), (50) 
for nonsingular v vortices 

^JL{0p) ^ IJLO + X"i« cos(6'p - Oh) + xa2v su\{Op -Oh), (51) 
and for for nonsingular w and uvw vortices 

fi{Op) = ^Io + aiu, cos{Op - Oh) + a2w sin{Op - Oh), (52) 

where we introduce the angle Oh characterizing the di- 
rection of perpendicular component of magnetic field 
H^ = H±{cos0h,smOh). 

Then from the Fig. ([3]) we can see that for the model 
vortex core (|48l49p the isoenergetic lines are indeed given 
by Eas. (|50l51l52p with /3i 7^ and /32 = for singular 
vortices; div,a2v 7^ for nonsingular v vortex; a2w 7^ 
and aiw = for nonsingular w and uvw vortices. In 
case if iJ^ = we have that all the coefficients are zero 
c^iv,o:2v,aiw,a2w = and the spectrum transformation 
should be of the next order in Hj^, i.e. proportional to (q- 
H_l)2 ~ H\cos{2{0p — Oh)), which is also demonstrated 
in FigIS] (see the plot for non-singular vortices; H^ = 0). 

In principle the numerical calculations described above 
allow to consider the spectrum transformation for the ar- 
bitrary values of Zeeman shift. However, the qualitative 



expression (|43)) have an advantage that they also take 
into account the self-consistent transformation of gap 
function in magnetic field A = A(H). If this transfor- 
mation is small it does not change the vortex symmetry 
therefore the general form of the energy spectrum con- 
serves yet with modified coefficients ai^2 in Eq.(|43|). In 
our further considerations we will take these coefficients 
as phcnomcnological constants since their particular val- 
ues does not affect the results qualitatively. 



and u vortices 



non-singular vortices; H =0. 

1i ' 1 



-1' ' 

0.5 

V vortex; H # 0. 



' 0.5 1 

w and uvw vortices; H # 0. 




0.5 

(ep-e^)/(27i) 




(ep-e,)/(23i) 



FIG. 3: Isoenergetic lines i-i{9p) corresponding to the differ- 
ent types of vortices in He3-B. Dash lines show fJ.{9p) in zero 
magnetic field and solid lines correspond to the modification 
of iJ.(6p) by the applied field. 



C. Magnetic susceptibility and energy dissipation. 

Now having in hand the expression for quasiparticle 
energy (^5)) it is then straightforward to calculate the 
non-stationary magnetization of vortex cores with the 
help if kinetic equation p6|) and the general expression 
for the magnetization (|32p . 

We will calculate the response magnetization of a single 
vortex using a unified expression for the spectrum (P5)) 
assuming that only the component of magnetic field per- 
pendicular to the vortex line varies in time. Also without 
loss of generality we put ai 7^ and a2 = 0, which can 
always be done by rotating the coordinate frame around 
the vortex axis by an appropriate angle. 

At first we will assume that an ac component of 
magnetic field is the one perpendicular to the vortex 
line. The component H, is assumed in general non- 
zero time independent. Also for the beginning the per- 
pendicular component is taken polarized along the x 
axis Hj^ ~ iJa;e*"'-'*x. Then solving the kinetic equa- 
tion (see Appendix B) and calculating the magnetization 
with the help of Eq. (jHS)) we obtain the following form of 
the magnetic susceptibility tensor Xxx — Xyy = X|| a-nd 
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Xxy 



-Xyx 



Xj_, where 



XII 



XX 



dkz a{ujrf + i/''')^rf 



dk. 



JyUJrf 



ixjy Lol - {urf + i/rY 



(53) 



(54) 



Here we have denoted LOy = —h~^den/diJ, which is the 
frequency corresponding to the interlevel energy spacing. 
Further wc will assume that this spacing is the same for 
all anomalous branches, i.e. cOy does not depend on the 
index n. Here wc have used the relation {hujy)dfo/ds ~ 
— 9/o/9/i and integrated over fi assuming for simplicity 
that LOy and a do not depend on fj. in order to perform the 
integration over /i. We have denoted a = N.i,a1/{8TT'^), 
where the overall factor Ny is a vortex density. 

The magnetic susceptibility defined by the Eqs. (j53l54p 
has resonances at the frequencies ujrf ~ ±cJt,. If the 
interlevel spacing cj„ depends on k^ [such as in Eq.®] 
the resonance peak transforms into the band with the 
absorption edge at w^/ = ojo = min|a;i,(fc2)|. In this 
case for the large enough relaxation times r 3> t^o the 
susceptibility is also peaked at w^/ = iwo behaving as 



X±,X|| 



1 



'rS 



The frequency independent magnitude of susceptibility 
can be estimated as 

XJ-,X|| ^ Xn^^v 

for singular vortices and 



X-L,X|| 



Xr 






e^. 



for nonsingular vortices, where Xn '^ i^oMB i^ ^ suscepti- 
bility of normal phase and A^ is the energy splitting of 
anomalous branches which can be taken of the order of 
bulk value of energy gap A. 

The resonant frequency of paramagnetic response 
where the absorption maximum takes place is selective 
to the polarization of external magnetic field. To demon- 
strate this let us consider two limiting cases: the lin- 
ear polarization H = i/j.e*'^'^^*x and the circular one 
H = i7_Le'"'-^*(x + iPy), where P = ±1 determines the 
direction of field rotation. In former case the dissipation 
determined by the Ea.(p3l) has the form 



Qe 



Urf 



H^ 



-Imx\\, 



where 



Imxw 



dk. 



aUJyUJrfT 



aU!y(jJrfl 



{Ujy 



UJrfYr'^ 



)2t2 + 1 (w, + CJ,/)2t2 + 1 



So the resonant energy absorption takes place both at 

LJrf = tdy and LJrf = —^y ■ 

In case of the circular polarization the dissipation rate 
is given by 



Qe 



Orf 



Hi 



Im{x\\ +ix±), 



where 



Im{x\\ +iX±) 



dk. 



UyUJj^fl 



— kf 



[uJy + PiJrfYT^ + I 



For the circular polarization the resonant frequency in 
the expression above depends on the direction of mag- 
netic field rotation w,. j = Pujy . 

Thus we see that the paramagnetic susceptibility as 
well as the dissipation losses in the system have the reso- 
nances defined in general by the interlevel spacing in the 
spectrum of vortex core fermion states (|41l42p . 

Interestingly, if we assume that the projection of mag- 
netic field onto the vortex line is also time-dependent 
Hz ~ H^ + Hz cos{uJrft + S) then besides the de- 
scribed above resonance at the main frequency ujrf — 
zth~^de/diJ, there can also appear resonances at higher 
frequencies w^/ = nujy and at fractional frequencies 
iJrf — idy/n, where n is an integer number. This sit- 
uation can naturally be realized in the experiment since 
it is not the real magnetic field which determines mag- 
netic response of vortices but a rotated one according to 
the Eq. (|T7)) . In general the rotation matrix R is spatially 
dependent hence the effective field HH has different di- 
rections at different points of the superfluid. To show 
the possibility of the additional resonances let us note 
that the expression for spectrum (j43p is in general non- 
linear in H. In particular, for the nonsingular vortices 
the coefficients in (|43p are proportional to H^ and there- 
fore depend on time as ai^2 = o;? 2 + 0^1,2 cos{u!rft + S). 
This additional modulation of the coefficients leads to the 
appearance of 2u!rf frequency terms in the energy spec- 
trum (|43|) . Then solving the kinetic equation in standard 
way described in Appendix ([B]) yields the resonances at 
Urf = ±ujy/2. In order to obtain the resonances at other 
frequencies it is necessary to consider the higher terms in 
spectrum expansion by the powers of the magnetic field 
H. 



IV. SUMMARY 

To sum up, we have investigated the spectrum of 
bound fermion states on vortices in "^He B modified by 
an external magnetic field. We have developed a general 
approach to study the spectrum perturbation based on 
symmetry grounds. It allowed us to determine qualita- 
tively up to the constants of the order unity the form of 
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bound states spectrum for different types of vortices in 
■^He B. An important advantage of this phenomenologi- 
cal approach to the spectral problem is a possibility to 
take into account the modification of the order parameter 
by the external field without extensive numerical calcu- 
lations. 

We also consider the paramagnetic susceptibility of 
fermionic ensemble bound within vortex cores. We have 
shown quite generally that it is the fermionic magneti- 
zation which determines the energy losses in ac external 
magnetic field driving the system out of equilibrium. It 
occurs that due to the coupling of orbital and spin quasi- 
particlc degrees of freedom the ac magnetic field induces 
transitions of bound fermions between different energy 
levels in a ladder of Caroli - de Gennes Matricon spec- 
trum ([1]). Consequently the paramagnetic susceptibility 
and energy absorption have resonance which occurs when 
the frequency of an ac external magnetic field equals the 
interlcvel spacing. Due to the broken time inversion sym- 
metry of vortex state the resonant behaviour of energy 
dissipation depends on the polarization of ac magnetic 
field. In particular for a circularly polarized magnetic 
field rotating over the vortex axis the presence of reso- 
nance depends on the relation between vortex winding 
direction and the direction of field rotation. 

Although the resonant absorption occurs at the same 
frequency for singular and nonsingular vortices the dis- 
sipation rate should be quite different in this two cases. 
Being proportional to Q^ ^ LOrf^^BH\ for the o and u 
vortices it is much less for the nonsingular vortices when 
Qs ~ Urf^J^BH\[nBH/ l^o)'^ since the magnetic field is 
assumed to much smaller than the spin depairing one so 
that /ifiiJ < Aq. 

Due to the resonant behaviour of paramagnetic suscep- 
tibility at the frequency of intcrlevcl transitions within 
the Caroli - de Gennes Matricon spectrum we can con- 
clude that measuring of a resonant magnetic susceptibil- 
ity of vortex cores in '^He B can provide a tool to study 
the discrete nature of bound fermions in vortex core. Also 
the difference in energy absorption rates for singular and 
non-singular vortices can provide an evidence for the par- 
ticular type of vortices realizing in "^He under different 
experimental conditions. 



Appendix A: Symmetries of quasiparticle spectrum. 



We are going to prove the general symmetries (|17I18|) 
of the spectrum of BdG system p3ll4p . At first let us 
consider the symmetry p7|) . Note that the coordinates 
in real space are related to the coordinates s, 6 as follows: 

X = s cos dry + b sin dn 



y 



9p — b cos 9p. 



Thus the transformation of s,b,Op to ~s,—b,6p + -k 
does not change the coordinates x^y. It means that 
the coordinate part of function Ak remains intact. At 
the same time changing kz,Op by — /c^j^p + tt leads to 
the reverse of momentum direction. Then the total 
gap function changes its sign under the transformation 
Ak(s, 6'p, fcz, 6) = — Ak(— s, 6'p + TT, — fc^, — 6). Note that 
the matrix A^ does not contain the Pauli matrix Uy 
therefore A^ 



Aj!. Then the complex conjugation of 



transformed BdG equations yields 



hk\ d 



AjcU* ^i-e-P]V 



m OS 



-e + P*\U 



(Al) 



(A2) 



Changing e by — e we obtain the system coinciding with 
the initial set of Eqs. (jl3ll4p which proves the relation 

m- 

Now let us consider the symmetry (|18p . We suppose 
here that Afc(H) = Aa;(-H). Note that the matrix A^ 
docs not contain the Pauli matrix ay therefore A^ — 
A J;. Then the complex conjugated BdG equations has 
the form 



m OS 



AkC/* = 



h^) 



V* 



i^^U* + k+V* = [e-P*]U* 
m OS 



(A3) 



(A4) 



Then changing H by — H we obtain initial set of 
Eas. (|13ll4p which proves (fT5|) . 
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Appendix B: Kinetic equation. 

In general the solution of kinetic equation (p9| can be 
found in the following form 



h^G{ep,t) 



de 



Then 



dt dt de de"^ 



(Bl) 
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where the second term can be neglected. Also, {H, G} = 
—ujyidG/dOp). Therefore for the functfon G{9p, t) we ob- 
tain the following equation: 



form: 



dG_ 

'dt 



dG G 



de 



(B2) 



where we denote hiUi, = —(de/dfi). 

Further we will consider the solution of kinetic equa- 
tion (|B2p when the quasiparticle spectrum is given by 
(|43| with ai 7^ 0, a2 = and Hj^ = H^e"^''f*x. Then we 
have 



iuJrfaiqxHxe''^'-'*-. 



de_ 

dt 



Let us search the solution of Eq. (|B2[) in the following 



G = (A cos 6'p + B sin 6lp)e" 



(B3) 



Substituting this form to the Eq. (|B2|) for the coefficients 
A and B we obtain: 



A ~ ujrfai 



B = iUrfCti 



(uJrf + i/T)H:c 



ujI - {uJrf + i/rf 



(B4) 



(B5) 



C. Caroli, P. G. de Gennes, J. Matricon, Phys. Lett. 9, 307 

(1964). 

G. E. Volovik, Pisma Zh. Exp. Teor. Fiz. 70, 601 (1999) 

[JETP Lett. 70, 609 (1999). 

N.B. Kopnin, Rep. Prog. Phys. 65, 1633 (2002). 

N.B. Kopnin and M.M. Salomaa, Phys. Rev. B 44, 9667 

(1991). 

G. E. Volovik, Zh. Eksp. Teor. Fiz. 104, 3070 (1993) 

[Sov. Phys. JETP 77, 435 (1993); G. E. Volovik, Pisma 

Zh. Eksp. Teor. Fiz. 57, 233 (1993) [JETP Lett. 57, 244 

(1993)]. 

M. Stone, Phys. Rev. B 54, 13222 (1996). 

N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

D.A. Ivanov, Phys. Rev. Lett. 86, 268 (2001); M. Stone, 

Suk-Bum Chung, Phys. Rev. B, 73, 014505 (2006). 

H. F. Hess, R. B. Robinson, R. C. Dynes, J. M. Valles, Jr., 

and J. V. Waszczak, Phys. Rev. Lett. 62, 214 (1989); H. 

F. Hess, R. B. Robinson, and J. V. Waszczak, Phys. Rev. 



Lett. 64, 2711 (1990). L Guillamon, H. Suderow, S. Vieira, 

L. Cario, P. Diener, and P. Rodiere, Phys. Rev. Lett. 101, 

166407 (2008). 
^° N.B. Kopnin and V.E. Kravtsov, Pisma Zh. Exp. Teor. Fiz. 

23, 631 (1976) [JETP Lett. 23, 578 (1976); N.B. Kopnin, 

7, 417 (1978) [JETP Lett. 7, 390 (1978). 
" N.B. Kopnin and G.E. Volovik, Phys. Rev. B 57, 8526 

(1997). 
^^ M.M. Salomaa and G.E. Volovik, Rev. Mod. Phys. 59, 533 

(1987). 
" G.E. Volovik and M.M. Salomaa, Phys. Rev. Lett. 51, 2040 

(1983). 
^* N.B. Kopnin Theory of nonequilibnum superconductivity 

Oxford University Press (2001). 
^^ R.M. White Quantum theory of magnetism Springer 

(1983). 
" M.A. Silaev, JETP Lett. 90, 433, (2009). 



